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^ ■ Abstract 

After a pedagogical introduction to the concept of resonance in classical and 
• quantum mechanics, some interesting applications are discussed. The subject in- 

P-h ! eludes resonances occurring as one of the effects of radiative reaction, the resonances 

■^j- ' involved in the refraction of electromagnetic waves by a medium with a complex re- 

fractive index, and quantum decaying systems described in terms of resonant states 
of the energy. Some useful mathematical approaches like the Fourier transform, the 
complex scaling method and the Darboux transformation are also reviewed. 



1 Introduction 
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Solutions of the Schrodinger equation associated to complex eigenvalues e = E — iY/2 
and satisfying purely outgoing conditions are known as Gamow-Siegert functions [1,2]. 
These solutions represent a special case of scattering states for which the 'capture' of 
the incident wave produces delays in the scattered wave. The 'time of capture' can be 
connected with the lifetime of a decaying system (resonance state) which is composed by 
the scatterer and the incident wave. Then, it is usual to take Re(e) as the binding energy 
of the composite while Im(e) corresponds to the inverse of its lifetime. The Gamow- 
Siegert functions are not admissible as physical solutions into the mathematical structure 
of Quantum Mechanics since, in contrast with conventional scattering wave-functions, 
they are not finite at r — > oo. Thus, such a kind of functions is acceptable in Quantum 
Mechanics only as a convenient model to solve scattering equations. However, because of 
the resonance states relevance, some approaches extend the formalism of quantum theory 
so that they can be defined in a precise form [3-8]. 

The concept of resonance arises from the study of oscillating systems in classical 
mechanics and extends its applications to physical theories like electromagnetism, optics, 
acoustics, and quantum mechanics, among others. In this context, resonance may be 
defined as the excitation of a system by matching the frequency of an applied force to a 
characteristic frequency of the system. Among the big quantity of examples of resonance 
in daily life one can include the motion of a child in a swing or the tuning of a radio 
or a television receiver. In the former case you must push the swing from time to time 
to maintain constant the amplitude of the oscillation. In case you want to increase the 
amplitude you should push 'with the natural motion' of the swing. That is, the acting of 
the force you are applying on the swing should be in 'resonance' with the swing motion. 
On the other hand, among the extremely large number of electromagnetic signals in space, 



your radio responds only to that one for which it is tuned. In other words, the set has to 
be in resonance with a specific electromagnetic wave to permit subsequent amplification 
to an audible level. In this paper we present some basics of resonance phenomenon. It is 
our intent to provide a strong primer introduction to the subject rather than a complete 
treatment. In the next sections we shall discuss classical models of vibrating systems 
giving rise to resonance states of the energy. Then we shall review some results arising 
from the Fourier transform widely used in optics and quantum mechanics. This material 
will be useful in the discussions on the effects of radiative reaction which are of great 
importance in the study of atomic systems. We leave for the second part of these notes 
the discussion on the resonances in quantum decaying systems and their similitudes with 
the behavior of optical devices including a complex refractive index. Then the complex 
scaling method arising in theories like physical chemistry is briefly reviewed to finish with 
a novel application of the ancient Darboux transformation in which the transformation 
function is a quantum resonant state of the energy. At the very end of the paper some 
lines are included as conclusions. 



2 Vibration, Waves and Resonances 
2.1 Mechanical Models 

Ideal vibrating (or oscillating) systems undergo the same motion over and over again. 
A very simple model consists of a mass m at the end of a spring which can slide back 
and forth without friction. The time taken to make a complete vibration is the period of 
oscillation while the frequency is the number of vibration cycles completed by the system 
in unit time. The motion is governed by the acceleration of the vibrating mass 

d 2 x ( k \ o /, \ 

-^ = -{m) X ^- W ° X (1) 

where wq ■— \Jkjm is the natural angular frequency of the system. In other words, a 
general displacement of the mass follows the rule 

x = A cos w t + B sin w t (2) 

with A and B two arbitrary constants. To simplify our analysis we shall consider a 
particular solution by taking A = a cos 6 and B = —a sin 9, therefore we can write 

x = a cos(w t + 6). (3) 

^ _ (2n+iyK-2e ^ n = Q ? 1 ; 2 . . ., the kinetic energy T = ^m(dx/ dt) 2 reaches its maximum 
value T max = mw\a 2 jl while x passes through zero. On the other hand, the kinetic 
energy is zero and the displacement of the mass is maximum [x = a is the amplitude 
of the oscillation) at t m = m ^~ e , m = 0, 1, 2, . . . This variation of T is just opposite of 
that of the potential energy V = kx 2 /2. As a consequence, the total stored energy E 
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is a constant of motion which is proportional to the square of the amplitude (twice the 
amplitude means an oscillation which has four times the energy): 

E = T + V = \mwla 2 . (4) 

Systems exhibiting such behavior are known as harmonic oscillators. There are plenty 
of examples: a weight on a spring, a pendulum with small swing, acoustical devices 
producing sound, the oscillations of charge flowing back and forth in an electrical circuit, 
the 'vibrations' of electrons in an atom producing light waves, the electrical and magnetical 
components of electromagnetic waves, and so on. 



2.1.1 Steady-state oscillations 

In actual vibrating systems there is some loss of energy due to friction forces. In other 
words, the amplitude of their oscillations is a decreasing function of time (the vibration 
damps down) and we say the system is damped. This situation occurs, for example, when 
the oscillator is immersed in a viscous medium like air, oil or water. In a first approach 
the friction force is proportional to the velocity Ff = — a4r, with a a damping constant 
expressed in units of mass times frequency. Hence, external energy must be supplied into 
the system to avoid the damping down of oscillations. In general, vibrations can be driven 
by a repetitive force F(t) acting on the oscillator. So long as F(t) is acting there is an 
amount of work done to maintain the stored energy (i.e., to keep constant the amplitude). 
Next we shall discuss the forced oscillator with damping for a natural frequency Wq and 
a damping constant a given. 

Let us consider an oscillating force defined as the real part of F(i) = Fe lwt = F e^ wt+V \ 
Our problem is to solve the equation 

d 2 x dx 2 „ [Fe iwt \ a 

~JPI + ^^7 + w o x = Re > 7:=—- 5 

dt 2 dt \ m J m 

Here the new damping constant 7 is expressed in units of frequency. The ansatz x = 
Ke(ze lwt ) reduces (|5]) to a factorizable expression of z, from which we get 

z = — , V = - )——. (6) 

m Wq — ttr + vyw 

We realize that z is proportional to the complex function Q, depending on the driving 
force's frequency w and parameterized by the natural frequency Wq and the damping 
constant 7. In polar form Q = |0|e^, the involved phase angle <fi is easily calculated by 
noticing that = e~ tl ^/\Q\ = — w 2 + i'yw, so we get 

tan0 = 9 7 9 . (7) 

Let us construct a single valued phase angle for finite values of Wq and 7. Notice that 
w < Wq leads to tan0 < while w — * Wq implies tan0 — > —00. Thereby we can set 
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(f)(w = 0) = and 4>(wq) = —n/2 to get G [— tt/2,0] for u> < wo- Now, since w > Wq 
produces tan0 > 0, we use tan(— </>) = — tan0 to extend the above defined domain 
G (— 7r, 0], no matter the value of the angular frequency w. Bearing these results in 
mind we calculate the real part of z (see equation ED to get the physical solution 



x = xq cos(wt + Tj + 



x :- 



F M 



m 



Notice that the mass oscillation is not in phase with the driving force but is shifted by 0. 
Moreover, 7^0 produces — > 0, so that this phase shift is a measure of the damping. 
Since the phase angle is always negative or zero (— n/2 < < 0), equation (JHJ) also means 
that the displacement x lags behind the force F(t) by an amount 0. On the other hand, 
the amplitude xq results from the quotient F /m scaled up by |0|. Thus, such a scale 
factor gives us a measure of the response of the oscillator to the action of the driving 
force. The total energy with a = Xq, is then a function of the angular frequency: 



E(w) 



(wqFq) 
2m 



\n\ 2 = 



(wqFq 
2m 



[Wn 



w 



2\2 



+ (7w) 2 _ 



(9) 



Equations ([7D and (jUJ comprise the complete solution to the problem. The last one, in 
particular, represents the spectral energy distribution of the forced oscillator with damping 
we are dealing with. It is useful, however, to simplify further under the assumption that 
7 << 1. For values of w closer to that of Wq the energy approaches its maximum value 
2mE(wo) ~ (Fo/7) 2 while E(w — > +00) goes to zero as w -4 . In other words, E(w) shows 
rapid variations only near wq. It is then reasonable to substitute 



W 



(wq - w)(w +w) « (wq - w)2w 



in the expressions of the energy and the phase shift to get 



E(w -> w ) w ( ~ ) ^O> w o,7), 



tan 1 



7 



2(w-w ) 



with 



(7/2) 2 



(10) 



'ID 



(12) 



(w -^) 2 + (7/2) 2 ' 

Equation ffl2|) describes a bell-shaped curve known as the Cauchy (mathematics), Lorentz 
(statistical physics) or Fock-Breit-Wigner (nuclear and particle physics) distribution. It 
is centered at w = Wo (the location parameter), with a half-width at half-maximum 
equal to 7/2 (the scale parameter) and amplitude (height) equal to 1. That is, the 
damping constant 7 defines the width of the spectral line between the half-maximum 
points Wo — w = ±7/2. Fig. [1] shows the behavior of the curve u for different values of 
the damping constant (spectral width) 7. 

These last results show that the adding of energy to the damping oscillator is most 
efficient if the vibrations are sustained at a frequency w = wq. In such a case it is said 
that the driving force is in resonance with the oscillator and wq is called the resonance 
frequency. Besides the above discussed spring-mass system, the motion of a child in a 
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wq — 1/2 wq wq + 1/2 



Figure 1: The Fock-Breit-Wigner (Lorentz-Cauchy) distribution u> for different values of the width at 
half maximum 7. 

swing is another simple example giving rise to the same profile. To keep the child+swing 
system oscillating at constant amplitude you must push it from time to time. To increase 
the amplitude you should push 'with the motion': the oscillator vibrates most strongly 
when the frequency of the driving force is equal to the frequency of the free vibration of 
the system. On the other hand, if you push against the motion, the oscillator do work 
on you and the vibration can be brought to a stop. The above cases include an external 
force steading the oscillations of the system. This is why equations ([THS]) are known as 
the steady- state solutions of the problem. 

We can calculate the amount of work W^ork which is done by the driving force. This 
can be measured in terms of the power P, which is the work done by the force per unit 
time: 

d TTT rl/ s dx dE ( dx\ 

P = _ H / TOk = F( ( )- = - + m7 y . (13) 

The average power (P) corresponds to the mean of P over many cycles. To calculate it we 
first notice that (dE/dt) = 0. That is, the energy E does not change over a period of time 
much larger than the period of oscillation. Now, since the square of any sinusoidal function 
has an average equal to 1/2, the last term in (fT3"|) has an average which is proportional to 
the square of the frequency times the amplitude of the oscillation. From ([H]) we get 

(P) = miW * Xl - (14) 

It is then clear that the driving force does a great deal of work to cancel the action of the 
friction force. In a similar form we obtain the average of the stored energy: 

(E) = mx o^ + w o\ (15 ) 

Remark that the mean of E does not depend on the friction but on the angular frequency 
of the driving force. If w is close to the resonance frequency wo, then (E) goes to the 
ideal oscillator's energy (j3J), scaled by (xo/a) 2 . Moreover, the same result is obtained no 
matter the magnitude of the driving force, since it does not play any role in (jT5l) . 
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2.1.2 Transient oscillations 

Suppose a situation in which the driving force is turned off at a given time t = to- This 
means no work is done to sustain the oscillations so that there is no supplied energy to 
preserve the motion any longer. This system can be studied by solving (jSJ) with F = 0. 
After introducing the ansatz x = Re(ze lwt ) we get a quadratic equation for w, the solution 
of which reads 

w± = i-f/2±<&, $ := sjw 2 - (7/2) 2 . (16) 

If 7 < wo then $ £ R and any of these two roots produces the desired solution: 

x = Re(ze iwt ) = \z\e-& cos{M + z ), t > t . (17) 

First, notice that the energy is not a constant of motion but decreases in exponential form 
E oc |,2| 2 e _7 *. The damping constant 7 is then a measure of the lifetime of the oscillation 
because at the time r = I/7, the energy is reduced to approximately the 36% (E — ► Eje) 
while the amplitude goes to the 60% of its initial value (\z\ — > \z\/y/e). Thus, the smaller 
the value of 7 the larger the lifetime r of the oscillation. In this way, for values of 7 
such that wo » 7/2, the discriminant in (fTBT) becomes 1? ~ wq. Thereby, the system 
exhibits an oscillation of frequency close to the resonance frequency wq. This means that 
large lifetimes are intimately connected with resonances for small values of the damping 
constant. 

As we can see, the resonance phenomenon is a characteristic of vibrating systems even 
in absence of forces steading the oscillations. Solutions like (fTTI) are known as transient 
oscillations because there is no force present which can ensure their prevalentness. They 
are useful to describe mechanical oscillators for which the driven force has been turned 
off at the time t = to or, more general, decaying systems like the electric field emitted 
by an atom. In general, 'resonance' is the tendency of a vibrating system to oscillate at 
maximum amplitude under certain frequencies w n , n — 0, 1, 2, ... At these resonance fre- 
quencies even small driving forces produce large amplitude vibrations. The phenomenon 
occurs in all type of oscillators, from mechanical and electromagnetic systems to quantum 
probability waves. A resonant oscillator can produce waves oscillating at specific frequen- 
cies. Even more, this can be used to pick out a specific frequency from an arbitrary 
vibration containing many frequencies. 

2.2 Fourier Optics Models 

In this section we shall review some interesting results arising from the Fourier transform. 
This mathematical algorithm is useful in studying the properties of optical devices, the 
effects of radiative reaction on the motion of charged particles and the energy spectra of 
quantum systems as well. Let {e lkx } be a set of plane waves orthonormalized as follows 

(e iKX , e ikx ) = lim f e l(k ~ K)x dx = lim 2 sm ^~ K)a = 2n6(k - k) (18) 
J _ a k — k 

with 5(x — xq) the Dirac's delta distribution. This 'function' arises in many fields of study 
and research as representing a sharp impulse applied at xq to the system one is dealing 
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with. The response of the system is then the subject of study and is known as the impulse 
response in electrical engineering, the spread function in optics or the Green's function in 
mathematical-physics. Among its other peculiar properties, the Dirac function is defined 
in such a way that it can sift out a single ordinate in the form 

/oo 
S(x-x )f(x )dx. (19) 
-oo 

In general, a one-dimensional function (p(x) can be expressed as the linear combination 

i r+oo 

(p{x) = — / ${k)e~ ikx dk (20) 



2vr _ 

where the coefficient of the expansion <p(k) is given by the following inner product 



1 

(e~ ikx , <p) = I ip(x)e ikx dx = — 
-oo 2tt 



+oo 



+oo 

( 

oo 



rOO 



<p(rv)dK 

(21) 



S(k — K)if(K)dK = <f(k). 



If (1201) is interpreted as the Fourier series of <p(x), then the continuous index k plays the 
role of an angular spatial frequency. The coefficient (p(k), in turn, is called the Fourier 
transform of <fi(x) and corresponds to the amplitude of the spatial frequency spectrum 
of <p(x) between k and k + dk. It is also remarkable that functions ff20|) and fT2"Tj) are 
connected via the Parseval formula 

\<p(x)\ 2 dx= / \p(k)\ 2 dk. (22) 

) J — oo 

This expression often represents a conservation principle. For instance, in quantum me- 
chanics it is a conservation of probability [9]. In optics, it represents the fact that all the 
light passing through a diffraction aperture eventually appears distributed throughout the 
diffraction pattern [10]. On the other hand, since x and k represent arbitrary (canonical 
conjugate) variables, if (p were a function of time rather than space we would replace x 
by t and then k by the angular temporal frequency w to get 

(p(t) = — / (p(w)e~ iwt dw, tp(w) = / if(t)e iwt dt (23) 

and 



DC 



\f{t)\ 2 dt= / \fi{w)\ 2 dw. (24) 

) J — oo 

Now, let us take a time depending wave <p(t), defined at x = by 



<p(t) = (poQ(t)e'^ cosw t, Q(t) = < (25) 
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Function (|25|) is a transient oscillation as it has been denned in the above sections. From 
our experience with the previous cases we know that it is profitable to represent (f(t) in 
terms of a complex function. In this case ip = Ke(Z), with 

Z(t) = A(t)e- iwot , A(t) = <p Q(t) exp (-It) . (26) 

Observe that |Z(t)| 2 = |A(t)| 2 . Then the Parseval formula (|24p gives 

+00 p+00 

\A(t)\ 2 dt= / \A(w)\ 2 dw. (27) 

00 J —00 

Since the stored energy at the time t is proportional to |A(t)| 2 , both integrals in equation 
( 1271) give the total energy W of the wave as it is propagating throughout x — 0. Thereby 
the power involved in the oscillation as a function of time is given by P(t) = dW/dt oc 
|v4(t)| 2 . In the same manner I w = dW/dw cx |A(w)| 2 is the energy per unit frequency 
interval. Now, from ([26]) we get A(t) = Z{t)e iw °\ so that 



(28) 



1 p+co ^ 1 r+00 

A(t) = — / Z(w)e- {w - Wo)t dw = — / a{e)e' i£t de 
2tt J^oo 2tt 

where e := w — w and a(e) := Z(e + w ). This last term is given by 

/+00 p+00 
A(t)e i£t dt = <p / e^-'^dt = rP r. (29) 
Jo i- i{w ~ w ) 



Then, up to a global constant, we have 



2(^o \ 2 (l/2f 



7 / (w - w ) 2 + (7/2) 2 ' 



(30) 



\<Po\ 



7=0.15 



7=0.3" 



1 y y ""^ — - vy -- 

t — » 

Fi gure 2: The power P(i) c>c |yo| 2 e 7 * involved with the transient oscillation |25|) for the same values 
of 7 as those given in Figure [TJ Remark that the area under the dotted curves is larger for smaller values 
of the line breadth 7. 

That is, the Fock-Breit-Wigner (FWB) function (|30|) is the spectral energy distribution 
of the transient oscillation (J25l) . In the previous section we learned that the inverse of 
the damping constant 7 measures the oscillation lifetime. The same rule holds for the 
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power |A(t)| 2 = |v?o| 2 e 7 *, as it is shown in Fig. [2j Let us investigate the extremal case 
of infinite lifetimes. Thus we calculate I w in the limit 7^0. The result reads 

I w ^2n(^]5(w-w ), 7^0 (31) 



where we have used 



7 



5{x) = -\\m[——\. (32) 



7T e-»0 V X 2 + 6 2 



Infinite lifetimes (I/7 = r — > +00) correspond to spectral energy distributions of in- 
finitesimal width (7 — > 0) and very high height (^0/7)- As a consequence, the transient 
oscillation (125]) has a definite frequency (w — * wo) along the time. The same conclusion is 
obtained for steady-state oscillations by taking F 2 = 8my?o in equations (JTH - fl2"]) . In the 
next sections we shall see that this lifetime<->width relationship links bound and decaying 
energy states in quantum mechanics. 

As a very simple application of the above results let the plane wave (125]) be the 
electric field emitted by an atom. Then W corresponds to the total energy radiated per 
unit area perpendicular to the direction of propagation. Equation (1311 in turn, relates 
in quantitative way the behavior of the power radiated as a function of the time to 
the frequency spectrum of the energy radiated. To give a more involved example let 
us consider a nonrelativistic charged particle of mass m e and charge q e , acted on by an 
external force F. The particle emits radiation since it is accelerated. To account for this 
radiative energy loss and its effect on the motion of the particle it is necessary to add a 
radiative reaction force F rad in the equation of motion to get 

P + P ^P + (0jP =me P ,33) 

with c the speed of light. This last expression is known as the A braham-Lorentz equation 
and is useful only in the domain where the reactive term is a small correction, since the 
third order derivative term does not fulfill the requirements for a dynamical equation (see 
e.g., reference [11], Ch 16). Bearing this condition in mind let us investigate the effect of 
an external force of the form F = —m e w\r. The Abraham-Lorentz equation is written 

2g e 2 W 3 ^ a 12 _ 2 ^ 

' r = —r + w%r. (34) 



3m R c 3 / dt 3 dt 2 

For small values of the third order term one has -C-f ~ — u> 2 r. That is, the particle 
oscillates like a mass at the end of a spring with frequency w . Hence -p-f ~ —w^f, so 
that the problem is reduced to the transient equation 

d 2 - d _ ~ 2<7 2 u> 2 . . 

_ r + 7 _ r+rao V = 0, 7 = s y, (35) 

the solution of which has the form (I25H26I) with ip(t) replaced by r(t) and (fo by a constant 
vector r . To get an idea of the order of our approach let us evaluate the quotient j/wq. 
A simple calculation gives the following constant 

^ ^ n «0/i 1 n— 23 



W 2 



0.624 x lQ-^s. (36) 
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Thus the condition 7 << Is -1 defines the appropriate values of the frequency wq. For 
instance, let 7 take the value 10 _3 s~ 1 . The appropriate value of wq is then of the order of 
an infrared frequency wq ~ 10 10 s _1 . But if 7 ~ lO -1 ^ -1 then wo ~ 10 6 s -1 , that is, the 
particle will oscillate at a radio wave frequency. 

Under the limits of our approach the radiative reaction force _F ra d plays the role of a 
friction force which damps the oscillations of the electric field. The resonant line shape 
defined by the Fock-Breit-Wigner function (13"U|) is broadened and shifted in frequency 
due to the reactive effects of radiation. That is, because the decaying of the power 
radiated P(t) oc |</9o| 2 e~ 7 *, the emitted radiation corresponds to a pulse (wave train) with 
effective length A ~ c/7 and covering an interval of frequencies equal to 7 rather than 
being monochromatic. The infinitesimal finiteness of the width in the spectral energy 
distribution ( 13T1) is then justified by the 'radiation friction'. In the language of radiation 
the damping constant 7 is known as the line breadth. 

Finally, it is well known that the effects of radiative reaction are of great importance in 
the detailed behavior of atomic systems. It is then remarkable that the simple plausibility 
arguments discussed above led to the qualitative features derived from the formalism of 
quantum electrodynamics. By proceeding in a similar manner, it is also possible to verify 
that the scattering and absorption of radiation by an oscillator are also described in terms 
of FBW-like distributions appearing in the scattering cross section. The reader is invited 
to review the approach in classical references like [11,12]. 



2.3 Fock's Energy Distribution Model 

Let {4>e(x)} be a set of eigenfunctions belonging to energy eigenvalues E in the continuous 
spectrum of a given one- dimensional Hamiltonian H. The vectors are orthonormalized as 
follows 

roo 

4> E (x)(j)E>(x)dx = 5(E' - E). (37) 



Notice we have taken for granted that the continuous spectrum is not degenerated, oth- 
erwise equation (1371) requires some modifications. Let us assume that a wave function 
ipo(x) can be expanded in a series of these functions, that is: 



M%)= / C{E)(j> E {x)dE, C(E)= / <j> E (x)ifo(x)dx. (38) 
The inner product of ipo with itself leads to the Parseval relation 

r+oo r+oo 

W:=^oM= \Mx)\ 2 dx= \C(E)\ 2 dE. (39) 



In the previous sections we learned that W allows the definition of the energy distribution 
00(E). In this case we have 

u{E):= d ^ = \C{E)\\ (40) 
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At an arbitrary time t > the state of the system reads 



ijj t {x) 



C{E)<p E {x)e- iEt/h dE. 



(41) 



The transition amplitude T(t > 0) from the state ipo m to ipt is given by the inner product 



T(t>G) = (M t ) 



ip (x)ip t (x)dx 



u){E)e- iEt/h dE. 



(42) 



Function T rules the transition probability |T(t)| 2 from ipo{x) to ipt{x) by relating the 
wave function at two different times: to an d t > t Q . It is known as the propagator in 
quantum mechanics and can be identified as a (spatial) Green's function for the time- 
dependent Schrodinger equation (see next section). From (|42|) . it is clear that T can be 
investigated in terms of spatial coordinates x or as a function of the energy distribution. 
Next, following the Fock's arguments [9], we shall analyze the transition probability for 
a decaying system by assuming that u(E) is given. Let u(E) be the Fock-Breit-Wigner 
distribution 



a; (a) 

Then equation ( 1421) reads 
T(t > 0) : 



1 


r (r/2) 2 " 


1 


(r/2) 2 j 


71 


_a? + (r/2) 2 _ 


71 


_(a + iT/2)(a-iT/2)_ 



a :- 



E — En. 



(43) 



1 

7T 



-iEot/h 



-iat/h 



(a + ir/2)(a-ir/2)' 



(44) 



1 

-R 


Im(z) 
► ;T 

R Re(z) 


\ t 


t-iT J 



Figure 3: Contour of integration in the complex i?-plane. 

We start with the observation that function ( 14*31) has two isolated singularities at a = 
±ir/2 (none of them lies on the real axis!). When R > 0, the point — iT/2 lies in the 
interior of the semicircular region which is depicted in Fig. [31 At this stage, it is convenient 
to introduce the function 

( \ -izt/n 

= TTTfW s(z) = (45) 
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Integrating f(z) counterclockwise around the boundary of the semicircular region the 
Cauchy integral formula (see e.g., [13]) is written 



I f(a)ala+ [ f(z)dz = 2mg(-iT/2). (46) 
Jr Jc r 



This last expression is valid for all values of R > 0. Since the value of the integral on the 
right in (j4"6"]) tends to as R — ► +oo, we finally arrive at the desired result: 



L /(G)rfG = T exp ("^)- 



(47) 



In this way the transition rate T(t) acquires the form of a transient oscillation 



r 

T(t > 0) = - exp(-iet/h), e := E - zT/2. (48) 

It is common to find T(t) as free of the factor T/2 (see e.g., [9], pp 159). This factor 
arises here because f l4"3]) has been written to be consistent with the previous expressions 
of a FWB distribution. The factor is easily removed if we take T/2 rather than (r/2) 2 
in the numerator of u(a). Now, we know that transient oscillations involve lifetimes and 
the present case is not an exception. Equation (148j) means that the transition from ipo(x) 
to ipt{x) is an exponential decreasing function of the time. Since this rate of change is 
symmetrical, it also gives information about the rate of decaying of the initial wave. Thus, 
the probability that the system has not yet decayed at time t is given by 

r x2 



\T(t)\= (jj exp(-rt/ft) s t>0. (49) 

We note that the state of the undecayed system ip does not change but decays suddenly. 
That is, the time t in fj4"9]) is counted off starting from the latest instant when the system 
has not decayed. The above description was established by Fock in his famous book on 
quantum mechanics [9] and this is why functions like (14"3"1) bears his name. It is remarkable 
that the first Russian edition is dated on August 1931. 

Finally, remark we have written ( 1481) in terms of the complex number e = Eq — iT/2. 
The reason is not merely aesthetic because, up to a constant factor, T(t > 0) is the Fourier 
transfornQ of the expansion coefficient C(E): 

C(E) := + = e) (50) 

where we have used (j4"3l . The relevance of this result will be clarified in the sequel. 



1 The Fourier transform in this case corresponds to the equations ([23]) . with the energy E and the 
angular frequency w related by the Einstein's expression E = hw. 
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3 Quanta, Tunneling and Resonances 



Let us consider the motion of a particle of mass m constrained to move on the straight-line 
in a given potential U(x). Its time- dependent Schrddinger equation is 



Hi/>(x,t) :-- 



2m dx 2 



+ U{x) 



d 

i/>(x,t) = ih—ip(x,t). 



(51) 



Let us assume the wave-function ip(x,t) is separable, that is ip(x,t) = <p(x)9(t). A simple 
calculation leads to 9(f) = exp(—iEt/h), with E a constant, and ip(x) a function fulfilling 



2m dx 2 



+ U{x) 



(f(x) = Hip(x) = Eip(x) 



(52) 



This time-independent Schrddinger equation (plainly the Schrddinger equation) defines a 
set of eigenvalues E and eigenf unctions of the Hamiltonian operator H which, in turn, 
represents the observable of the energy. Now, to get some intuition about the separableness 
of the wave-function let us take the Fourier transform of its temporal term 



9(E) 



lim /" e i{E - E ' )t/h dt = 2it5(E - E') 

J —a 



(53) 



where we have used f ll8p . This last result means that the energy distribution is of in- 
finitesimal width and very high height. In other words, the system has a definite energy 
E = E' along the time. Systems exhibiting this kind of behavior are known as stationary 
and it is said that they are conservative. Since the Hamiltonian operator H does not 
depend on t and because for any analytic function / of H one has 



f(HMx) = f(E)(p(x) 
our separability ansatz ip — >• tp9 can be written 



ip(x, t) = exp 



(p(x) 



e~ iEt / h V (x) 



(54) 



(55) 



According with the Born's interpretation, the wave function <p(x) defines the probability 
density p(x) = \<p(x)\ 2 of finding the quantum particle between x and x + dx. Thereby, 
the sum of all probabilities (i.e., the probability of finding the particle anywhere in the 
straight-line at all) is unity: 



+oo 



p(x)dx 



ip(x)ip(x)dx 



\ip(x)\ 2 dx 



(56) 



The above equation represents the normalization condition fulfilled by the solutions <p(x) 
to be physically acceptable. Hence, they are elements of a vector space Ti, consisting of 
square-integrable functions and denoted as H = L 2 (H,fi), with \x the Lebesgue measure 
(for simplicity in notation we shall omit [i by writing L 2 (R)). As an example, TC can be the 
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space spanned by the Hermite polynomials H n (x), weighted by the factor fi(x) = e x I 2 
and defined as follows 

H n (x) = (-l)« e * 2 /2j_ e -* 2 /2. (57) 

In quantum mechanics, observables are represented by the so-called Hermitian operators 
in the Hilbert space 7i. A differential operator A defined on L 2 {R) is said to be Hermitian 
if, whenever Af and Ag are defined for f,g G L 2 (R) and belong to L 2 (R), then 

/+oo _ 
f(x)Ag(x)fi(x)dx. (58) 
-oo 

In particular, if / = g the above definition means that the action of the Hermitian operator 
A on / G L 2 (R) is symmetrical. If Af(x) = af(x), we have (Af, f) = a and (/, Af) = a, 
so that a = a =>■ a G R. In other words, the eigenvalues of a Hermitian operator acting 
on L 2 (H) are real numbers. It is important to stress, however, that this rule is not true 
in the opposite direction. In general, as we are going to see in the next sections, there is 
a wide family of operators A\ sharing the same set of real eigenvalues {a}. (The family 
includes some non- Hermitian operators!) Moreover, notice that the rule is not necessarily 
true if f(x) £ L 2 (R). In general, the set of solutions of Af(x) = af(x) is wider than 
L 2 (R). That is, the complete set of mathematical solutions <p(x) embraces functions such 
that its absolute value |y(x)| diverges even for real eigenvalues. A plain example is given 
by the solutions representing scattering states because they do not fulfill f|56l) . In such a 
case one introduces another kind of normalization like that defined in ( l3"Tj) . with a similar 
notion of the Born's probability as we have seen in the previous section. Normalization 
( I56p is then a very restrictive condition picking out the appropriate physical solutions 
among the mathematical ones. This is why the Schrodinger equation (j5ip is "physically 
solvable" for a very narrow set of potentials. 

If E is real, the time-dependent factor in (155]) is purely oscillatory (a phase) and 
the time displacement ip(x,t) gives the same 'prediction' (probability density) as <p(x). 
As a result, both of these vectors lead to the same expectation values of the involved 
observables 

/+oo f+OO 
ip(x,t)Aip(x,t)dx = / lp(x)A(p(x) = a. (59) 
-oo J — oo 

In particular, (H) = E shows that the eigenvalue E is also the expectation value of 
the energy. Notice that stationary states are states of well-defined energy, E being the 
definite value of its energy and not only its expectation value (see equation [53]). That 
is, any determination of the energy of the particle always yields the particular value E. 
Again, as an example, let us consider a scattering state. Far away from the influence of 
the scatterer, it is represented by a plane wave like 

i/j(x,t) = exp ^— — Et^j exp ^— j^V^j (60) 

where p is the linear momentum of the particle. Function (160]) represents a state having a 
definite value for its energy. However, there is no certainty neither on the position of the 
particle nor in the transit time of the particle at a given position. In general, the energy 
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distribution will not be a continuous function. It could include a set of isolated points 
(discrete energy levels) and/or continuous portions showing a set of very narrow and high 
peaks (resonance levels). The former correspond to infinite lifetime states (the observed 
discrete energy levels of atoms are good examples) while the lifetime of the second ones 
will depend on the involved interactions. 



3.1 Quasi-stationary States and Optical Potentials 



It is also possible to define a probability current density j 

h 



J 



2mi 



-( d 

^Vdx 4 



(61) 



which, together with the probability density p = \ip\ 2 satisfies a continuity equation 

dt dx 



(62) 



exactly as in the case of conservation of charge in electrodynamics. Observe that station- 
ary states fulfill dp/dt = 0, so that p ^ p(t) and j = 0. What about decaying systems for 
which the transition amplitude T(t > 0) involves a complex number e = Eq — iT/2 like 
that found in equation (HH!) ? Let us assume a complex eigenvalue of the energy Hip € = e<p e 
is admissible in (1551) . Then we have p(x,t) = p e (x)e~ rt ^ h and j ^ 0. That is, complex 
energies are included at the cost of adding a non-trivial value of the probability current 
density j. A conventional way to solve this 'problem' is to consider a complex potential 
U = Ur + iUi. Then equation (1521 acquires the form 



(63) 



the integration of which can be identified with the variation of the number of particles 



dt h 



Ui(x)p(x)dx. 



(64) 



Let Ui(x) = Uq be a constant. If Uq > 0, there is an increment of the number of particles 
(dN/dt > 0) and viceversa, Uo < leads to a decreasing number of particles (dN/dt < 0). 
In the former case the imaginary part of the potential works like a source of particles 
while the second one shows Uj as a sink of particles. The introduction of this potential 
into the Schrodinger equation gives 



T?_cP_ 
2m dx 2 



+ U R (x) + iU 



<p e (x) 



En 



ip e [x) 



(65) 



Then, the identification Uq = —T/2 reduces the solving of this last equation to a stationary 
problem. The exponential decreasing probability density p e {x)e~ vt ^ h is then justified by 
the presence of a sink-like potential Ui(x) = Uo < 0. However, this solution requires the 
introduction of a non-Hermitian Hamiltonian because the involved potential is complex. 
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Although such a Hamiltonian is not an observable in the sense defined in the above section, 
notice that U = —iT/2 is a kind of damping constant. The lifetime of the probability 
p(x, t) is defined by the inverse of Y and, according with the derivations of the previous 
section, the energy distribution shows a bell-shaped peak at E — E Q . In other words, the 
complex eigenvalue e = E — iT/2 is a pole of uu(E) and represents a resonance of the 
system. The modeling of decay by the use of complex potentials with constant imaginary 
part is know as the optical model in nuclear physics. The reason for such a name is 
clarified in the next section. 

3.1.1 Complex Refractive index in Optics 

It is well known that the spectrum of electromagnetic energy includes radio waves, infrared 
radiation, the visible spectrum of colors red through violet, ultraviolet radiation, x-rays 
and gamma radiation. All of them are different forms of light and are usually described as 
electromagnetic waves. The physical theory treating the propagation of light is due mainly 
to the work of James Clerk Maxwell (1831-1879). The interaction of light and matter or 
the absorption and emission of light, on the other hand, is described by the quantum 
theory. Hence, a consistent theoretical explanation of all optical phenomena is furnished 
jointly by Maxwell's electromagnetic theory and the quantum theory. In particular, the 
speed of light c = 299, 792,456.2 ± 1.1m/ s is a part of the wave equation fulfilled by the 

— * — * 

electric field E and the magnetic field H: 

A = E,H. (66) 

The above expression arises from the Maxwell's equations in empty space with c = 
(yUoeo) -1 / 2 ■ The constant /io is known as the permeability of the vacuum and the con- 
stant Cq is called the permitivity of the vacuum [11]. In isotropic nonconducting media 
these constants are replaced by the corresponding constants for the medium, namely 
p, and e*. Consequently, the speed of propagation v of the electromagnetic fields in a 
medium is given by v = (/ic*) -1 / 2 . The index of refraction n is defined as the ratio of 
the speed of light in vacuum to its speed in the medium: n — c/v. Most transparent 
optical media are nonmagnetic so that /i//io = 1, in which case the index of refraction 
should be equal to the square root of the relative permitivity n = (e*/eg) 1 / 2 = K x l 2 . In 
a nonconducting, isotropic medium, the electrons are permanently bound to the atoms 
comprising the medium and there is no preferential direction. This is what is meant by 
a simple isotropic dielectric such a glass [14]. Now, consider a (one-dimensional) plane 
harmonic wave incident upon a plane boundary separating two different optical media. In 
agreement with the phenomena of reflection and refraction of light ruled by the Huygen's 
principle, there will be a reflected wave and a transmitted wave (see e.g. [10]). Let the 
first medium be the empty space and the second one having a complex index of refraction 

M = n + ini. (67) 

Then, the wavenumber of the refracted wave is complex: 

IC = k + ia (68) 
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and we have 



(incident wave) 
(reflected wave) 
£n e -ax e i(kx-wt) ( re f rac t e d wave) 

If a > the factor e~ ax indicates that the amplitude of the wave decreases exponentially 
with the distance. That is, the energy of the wave is absorbed by the medium and varies 
with distance as e~ 2ax . Hence 2a is the coefficient of absorption of the medium. The 
imaginary part ni of Af, in turn, is known as the extinction index. In general, it can 
be shown that the corresponding polarization is similar to the amplitude formula for a 
driven harmonic oscillator [14]. Thus, an optical resonance phenomenon will occur for 
light frequencies in the neighborhood of the resonance frequency Wq = (K/m) 1 ^ 2 . The 
relevant aspect of these results is that a complex index of refraction M leads to a complex 
wavenumber /C. That is, the properties of the medium (in this case an absorbing medium) 
induce an specific behavior of the electromagnetic waves (in this case, the exponentially 
decreasing of the amplitude). This is the reason why the complex potential discussed in 
the previous section is named the 'optical potential'. 

3.2 Quantum Tunneling and Resonances 

In quantum mechanics the complex energies were studied for the first time in a paper 
by Gamow concerning the alpha decay (1928) [1]. In a simple picture, a given nucleus 
is composed in part by alpha particles (f-ffe nuclei) which interact with the rest of the 
nucleus via an attractive well (obeying the presence of nuclear forces) plus a potential 
barrier (due, in part, to repulsive electrostatic forces). The former interaction constrains 
the particles to be bounded while the second holds them inside the nucleus. The alpha 
particles have a small (non-zero) probability of tunneling to the other side of the barrier 
instead of remaining confined to the interior of the well. Outside the potential region, 
they have a finite lifetime. Thus, alpha particles in a nucleus should be represented by 
quasi-stationary states. For such states, if at time t = the probability of finding the 
particle inside the well is unity, in subsequent moments the probability will be a slowly 
decreasing function of time (see e.g. Sections 7 and 8 of reference [9]). In his paper of 
1928, Gamow studied the escape of alpha particles from the nucleus via the tunnel effect. 
In order to describe eigenf unctions with exponentially decaying time evolution, Gamow 
introduced energy eigenf unctions ipc belonging to complex eigenvalues Z G = E G — iT G , 
Tq > 0. The real part of the eigenvalue was identified with the energy of the system and 
the imaginary part was associated with the inverse of the lifetime. Such 'decaying states' 
were the first application of quantum theory to nuclear physics. 

Three years later, in 1931, Fock showed that the law of decay of a quasi-stationary 
state depends only on the energy distribution function u(E) which, in turn, is meromor- 
phic [9]. According to Fock, the analytical expression of u(E) is rather simple and has 
only two poles E = E ±iT, T>0 (see our equation (j4"3"|) and equation (8.13) of [9]). A 
close result was derived by Breit and Wigner in 1936. They studied the cross section of 



E = Eq e^Ooz-Mrf) 
E' = Eq Q i ( k 'o x ~ wt ) 
E" = En e*^' - ™*) = 
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slow neutrons and found that the related energy distribution reaches its maximum at Er 
with a half-maximum width Tr. A resonance is supposed to take place at Er and to have 
"half-value breath" Tr [2]. It was in 1939 that Siegert introduced the concept of a purely 
outgoing wave belonging to the complex eigenvalue e = E — iT/2 as an appropriate tool in 
the studying of resonances [15]. This complex eigenvalue also corresponds to a first-order 
pole of the S matrix [16] (for more details see e.g. [17]). However, as the Hamiltonian is 
a Hermitian operator, then (in the Hilbert space TC) there can be no eigenstate having a 
strict complex exponential dependence on time. In other words, decaying states are an ap- 
proximation within the conventional quantum mechanics framework. This fact is usually 
taken to motivate the study of the rigged (equipped) Hilbert space Ti [3-5] (For a recent 
review see [6]). The mathematical structure of 7i lies on the nuclear spectral theorem 
introduced by Dirac in a heuristic form [18] and studied in formal rigor by Maurin [19] 
and Gelfand and Vilenkin [20]. 

In general, solutions of the Schrodinger equation associated to complex eigenvalues 
and fulfilling purely outgoing conditions are known as Gamow- Siegert functions. If u is a 
function solving Hu = eu, the appropriate boundary condition may be written 

lim (u'^iku)= lim {{—(3 =F ik) u} = 0, (69) 

x— >±oo x— >±oo 

with (3 defined as the derivative of the logarithm of u: 

(3:=-^-\nu. (70) 

ax 

Now, let us consider a one-dimensional short-range potential U(x), characterized by a 
cutoff parameter ( > 0. The general solution of (l52~l) can be written in terms of ingoing 
and outgoing waves: 

U< := u(x < -C) = Ie ikx + Le~ ikx , U> := u(x >£) = Ne ~ ikX + Se ikx ( 71 ) 

where the coefficients /, L, N, S, depend on the potential parameters and the incoming 
energy k 2 = 2me/h 2 (the kinetic parameter k is in general a complex number k = kn+iki), 
they are usually fixed by imposing the continuity conditions for u and du/dx at the points 
x = ±£. Among these solutions, we are interested in those which are purely outgoing 
waves. Thus, the second term in each of the functions (|71|) must dominate over the first 
one. For such states, equation (IBTI) takes the form: 

j < = -t,| M< | 2 , j >=v \ u> \ 2 , v:= ^^ + k) = ^-. (72) 

2m m 

This last equation introduces the flux velocity v. If e is a real number e = E, then k 
is either pure imaginary or real according to E negative or positive. If we assume that 
the potential admits negative energies, we get k± = ±i^\2mE/h 2 \ and f!72|) vanishes 
(the flux velocity v oc k^ is zero outside the interaction zone). Notice that the solutions 
connected with k + , are bounded so that they are in L 2 (R). That is, they are the 
physical solutions associated with a discrete set of eigenvalues 2mE n /h 2 = k 2 + solving 
the continuity equations for u and du/dx at x — ±£. On the other hand, antibound states 
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u* E increase exponentially as \x\ — > +00. To exhaust the cases of a real eigenvalue e, let 
us take now 2mE/h 2 = k 2 > 0. The outgoing condition (1691) drops the interference term 
in the density 

p(x;t) = \N\ 2 + \S\ 2 + 2\NS\cos(2Kx + AigS/N) } x > (, (73) 

so that the integral of p = \S\ 2 is not finite neither in space nor in time (similar expressions 
hold for x < — £). Remark that flux velocity is not zero outside the interaction zone. 
Thereby, E > provides outgoing waves at the cost of a net outflow j ^ 0. To get 
solutions which are more appropriate for this nontrivial j, we shall consider complex 
eigenvalues e. Let us write 

e = E- l -T, e R <xk? R -k%, e I (x2k R k I (74) 

where 2me/h 2 = (k R + ikj) 2 . According to (1691) . the boundary condition for f3 reads now 

lim {-p ± (fcj - ik R )} = (75) 

x— >±oo 

so that the flux velocity is v + oc for x > ( and t>_ oc —k R for x < — £. Hence, the 
"correct" direction in which the outgoing waves move is given by k R > 0. In this case, 
the density 

p(z; t) = \u(x, t)\ 2 = e~ rt/h \u(x)\ 2 , lim p(x; t) oc e - r (*-^±)/^ ( 76 ) 

x— >±oo 

can be damped by taking T > 0. Thereby, kj ^ and k R ^ have opposite signs. Since 
/cr > has been previously fixed, we have kj < 0. Then, purely outgoing, exponentially 
increasing functions (resonant states) are defined by points in the fourth quadrant of the 
complex fc-plane. In general, it can be shown that the transmission amplitude S in (|71|) is 
a meromorphic function of k, with poles restricted to the positive imaginary axis (bound 
states) and the lower half-plane (resonances) [21]. Let k n be a pole of S in the fourth 
quadrant of the fc-plane, then — k n is also a pole while k n and — k n are zeros of S (see 
Figure HI left). On the other hand, if S is studied as a function of e, a Riemann surface 
of e 1//2 = k is obtained by replacing the e-plane with a surface made up of two sheets 
R and Ri, each cut along the positive real axis and with Ri placed in front of R (see 
e.g. [13], pp 337). As the point e starts from the upper edge of the slit in R and describes 
a continuous circuit around the origin in the counterclockwise direction (Figure HI right), 
the angle increases from to 2it. The point then passes from the sheet i?o to the sheet 
Ri, where the angle increases from 2ir to Air. The point then passes back to the sheet Ro 
and so on. Complex poles of S(e) always arise in conjugate pairs (corresponding to k and 
—k) while poles on the negative real axis correspond to either bound or antibound states. 

Observe that density ( 1761) increases exponentially for either large \x\ or large negative 
values of t. The usual interpretation is that the compound (ip e , V) represents a decaying 
system which emitted waves in the remote past t — x/v. As it is well known, the long 
lifetime limit (r — > 0) is useful to avoid some of the complications connected with the 
limit t — > —00 (see discussions on time asymmetry in [22]). In this context, one usually 
imposes the condition: 

51 « L < 77 > 
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Figure 4: Left: Schematic representation of the poles (disks) and the zeros (circles) of the transmission 
amplitude S(k) in the complex fc-plane. Bounded energies correspond to poles located on the positive 
imaginary axis Right: The two-sheet Riemann surface = k. The lower edge of the slit in Rq is joined 
to the upper edge of the slit in R±, and the lower edge of the slit in R\ is joined to the upper edge of the 
slit in i?o- The picture is based on the description given by J.M. Brown and R.V. Churchill in ref. [13]. 

Thus, the level width V must be much smaller than the level spacing AE in such a way that 
closer resonances imply narrower widths (longer lifetimes). In general, the main difficulty 
is precisely to find the adequate E and T. However, for one-dimensional stationary short 
range potentials, in [21] it has been shown that the superposition of a denumerable set 
of FBW distributions (each one centered at each resonance E n ,n = 1,2, . . .) entails an 
approximation of the coefficient T such that the larger the number N of close resonances 
involved, the higher the precision of the approximation (see Fig. [Sj): 

N 

T^uj N {e R )=^2uj{e R ,E n ) (78) 

n=l 

with 




10 20 30 40 50 



Figure 5: Functions T and u>n (dotted curve) for a square well with strength Vo = 992.25 and weigh 
b = 20 for which the FWB sum matches well the transmission coefficient for the first five resonances [23] 
(see also [21]). 

Processes in which the incident wave falls upon a single scatterer are fundamental in the 
study of more involved interactions [24]. In general, for a single target the scattering 
amplitude is a function of two variables (e.g. energy and angular momentum). The above 
model corresponds to the situation in which one of the variables is held fixed (namely, the 



20 



angular momentum). A more realistic three dimensional model is easily obtained from 
these results: even functions are dropped while an infinitely extended, impenetrable wall 
is added at the negative part of the straight line [25]. Such a situation corresponds to 
s- waves interacting with a single, spherically symmetric, square scatterer (see e.g. [26]). 

3.3 Complex Scaling Method 

Some other approaches extend the framework of quantum theory so that quasi-stationary 
states can be defined in a precise form. For example, the complex-scaling method [7,8,27] 
(see also [28]) embraces the transformation H — > UHU^ 1 = H e , where U is the complex- 
scaling operator U = e - dxp / h ) with 9 a dimensionless parameter and [X, P] = ih. The 
transformation is achievable by using the Baker-Campbell-Hausdorff formulae [29]: 



oo 1 

Be- A = {e A ,B} = J2-{A n ,B} (80) 



e A 

n=0 n! 



with A and B two arbitrary linear operators and 

{A n ,B} = [A,[A,...[A, B}...}}. 
n times 

The identification A = —9XP/ h and B = X leads to 

UXU- 1 = -{iO) n X = e w X, UPU- 1 = ^ -{-i9) n P = e^P, (81) 

n=0 n ' n=0 n ' 

where we have used 

{{XP) n , X} = {-ih) n X, {(XP) n , P} = (ih) n P. (82) 
The following calculations are now easy 

UP 2 U^ = (UPU' 1 ) (UPU- 1 ) = (UPU- 1 ) 2 = e- 2l9 P 2 , 



UV(X)U- 1 = U ( jr ljVkX k J U- 1 = V (e w X) 
\k=o ' / 



(83) 



So that we finally get 

UHU- 1 =H e = e - 2ie P 2 + V(e ie X). (84) 
Remark that in the Schrodinger's representation we have 

X = x, P = -ih-^, U = e ldxd/dx Uf(x) = f(xe td ). (85) 

This transformation converts the description of resonances by non-integrable Gamow- 
Siegert functions into one by square integrable functions. Let k = \k\e ia = \k\e~ llS be 
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lm(k) 




Re(k) 



Figure 6: Polar form of an arbitrary point on the complex /c-plane. 



a point on the complex fc-plane (see Figure [6]). If k lies on the fourth quadrant then 
< (3 < n/2 and the related Gamow-Siegert function u t behaves as follows 

U e (x -> ±00) ~ e ±<[fc|*«w/3 e ±|fc|x 8 ln/9 i (86) 

That is, u e diverges for large values of \x\. The behavior of the complex-scaled function 
u t = U(u e ), on the other hand, reads 

U € (X -> ±00) ~ e ±i\k\xcos(e-f3) e ^\k\x S in(0-f})_ ( g? ) 

Thereby, u e is a bounded function if 9 — (3 > 0, i.e., if tan# > tan (3. The direct calculation 
shows that complex-scaling preserves the square-integrability of the bounded states <p n , 
whenever < 9 < tt/2. Then one obtains 

< 9 - (3 < tt/2. (88) 




Fi gure 7: Left. Schematic representation of the two-sheet Riemann surface showed in FigureJH Dashed 
curves and empty squares lie on the first sheet Rq, continuous curves and fulled squares lie on the second 
sheet i?i Right. The complex rotated plane 'exposing' the resonant poles lying on the first sheet Rq. 



As regards the complex-scaled scattering states we have 

g±i[fc|a; ^ g±i|fe[a; cos 6^f\k\x sin£> (89) 

So that plane waves are transformed into exponential decreasing or increasing functions for 
large values of \x\. To preserve the oscillating form of scattering wave-functions the kinetic 
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parameter k has to be modified. That is, the transformation k = \k\ — > | A; | e — z6! reduces 
( 189]) to the conventional plane-wave form of the scattering states. This transformation, 
however, induces a rotation of the positive real axis in the clockwise direction by the angle 
28: E oc k 2 — > \k\ 2 e~ l2e oc Ee~ l2e . That is, the rotated energy is complex e = Er — iY/2 
with Er = Ecos(26), T/2 = sin(2#). In summary, complex rotation is such that: 1) 
Bound state poles remain unchanged under the transformation 2) Cuts are now rotated 
downward making an angle of 29 with the real axis 3) Resonant poles are 'exposed' by 
the cuts (see Figure [7]). Another relevant aspect of the method is that it is possible 
to construct a resolution to the identity [30]. Moreover, as the complex eigenvalues are 
^-independent, the resonance phenomenon is just associated with the discrete part of 
the complex-scaled Hamiltonian [31] (but see [28]). As a final remark, let us enphasize 
that complex-scaling 'regularizes' the divergent Gamow-Siegert functions u € at the cost 
of introducing a non-Hermitian Hamiltonian Hg. From equation we get 

Hi = e i2e P 2 + V(e~ l9 R) ^ H e . (90) 

In other words, the 'regularized' solutions u e are square-integrable eigenfunctions of a 
complex potential V(e td x) belonging to the complex eigenvalue e. 



3.4 Darboux-Gamow Transformations 

In a different survey, complex eigenvalues of Hermitian Hamiltonians have been used to 
implement Darboux (supersymmetric) transformations in quantum mechanics [32-37] (see 
also the discussion on 'atypical models' in [38]). The transformed Hamiltonians include 
non-Hermitian ones, for which the point spectrum sometimes has a single complex eigen- 
value [21,25,33-35]. This last result, combined with appropriate squeezing operators [39], 
could be in connection with the complex-scaling technique. In general, supersymmetric 
transformations constitute a powerful tool in quantum mechanics [38]. However, as far 
as we know, until the recent results reported in [21, 25, 35, 36] the connection between 
supersymmetric transformations and resonant states has been missing. In this context 
and to throw further light on the complex function (3 we may note that flTUj) transforms 
the Schrodinger equation (1521) into a Riccati one 

-(3' + (3 2 + e = V (91) 

where we have omitted the units. Remark that (I9TI) is not invariant under a change in 
the sign of the function (3: 

P> + pi + e = V + 2p'. (92) 

These last equations define a Darboux transformation V = V(x, e) = V(x) + 2(3' (x) of 
the initial potential V. This transformation necessarily produces a complex function if u 
in equation fl701 is a Gamow-Siegert function u e . That is, a Darboux-Gamow deformation 
is defined as follows [21]: 

d 2 

V = V + 2(3' = V -2— \nu e . (93) 
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The main point here is that the purely outgoing condition (|69|) leads to (3' — > so that 
V — > V, in the limit |x| — > +00. In general, according with the excitation level of the 
transformation function u e , the real Vr and imaginary Vi parts of V show a series of 
maxima and minima. Thus, the new potential behaves as an optical device emitting 
and absorbing probability flux at the same time, since the function Ii(x) shows multiple 
changes of sign [21,25,35,36]. On the other hand, the solutions y = y(x,e,£) of the 
non-Hermitian Schrodinger equation 

- y" + Vy = £y (94) 

are easily obtained 

V =c (95) 

where W(*, *) stands for the Wronskian of the involved functions and ip is eigen-solution 
of (1941) with eigenvalue £. It is easy to show that scattering waves and their Darboux- 
Gamow deformations share similar transmission probabilities [21]. Now, let us suppose 
that Hamiltonian H includes a point spectrum <Jd{H) C Sp(H). If ip n is a (square- 
integrable) eigenfunction with eigenvalue £ n , then its Darboux-Gamow deformation (1951) 
is bounded: 



lim y n = =F(y/£ n + ik)( lim ip n ). (96) 

x— >±oo x— >±oo 

Thereby, y n is a normalizable eigenfunction of H with eigenvalue £ n . However, as e is 
complex, although the new functions {y n } rnay be normalizable, they will not form an 
orthogonal set [34] (see also [40] and the 'puzzles' with self orthogonal states [41]). There 
is still another bounded solution to be considered. Function y e oc ip~ l fulfills equation (1941) 
for the complex eigenvalue e. Since lim^-too \y e \ 2 = e ±2klX and kj < 0, we have another 
normalizable function to be added to the set {y„}. 

In summary, one is able to construct non-Hermitian Hamiltonians H for which the 
point spectrum is also ad(H), extended by a single complex eigenvalue Sp(if) = Sp(-ff) U 
{e}. As we can see, the results of the Darboux-Gamow deformations are quite similar 
to those obtained by means of the complex-scaling method. This relationship deserves a 
detailed discussion which will be given elsewhere. 



4 Conclusions 

We have studied the concept of resonance as it is understood in classical mechanics by 
analyzing the motion of a forced oscillator with damping. The resonance phenomenon 
occurs for steady state oscillations when the driving force oscillates at an angular frequency 
equal to the natural frequency of the oscillator. Then the amplitude of the oscillation 
is maximum and wq is called the resonance frequency. The spectral energy distribution 
corresponds to a Fock-Breit-Wigner (FBW) function, centered at wq and having a line 
breadth equal to the damping constant 7. Resonance is present even in the absence of 
external forces (transient oscillations). In such case the energy decreases exponentially 
with the time so that the damping constant 7 is a measure of the lifetime of the oscillation 
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r = I/7. Similar phenomena occur for the electromagnetic radiation. In this context, the 
effects of radiative radiation can be approximated by considering the radiative reaction 
force friction force which damps the oscillations of the electric field. Thus, the 

concept of resonance studied in classical mechanics is easily extended to the Maxwell's 
electromagnetic theory. The model also applies in vibrating elastic bodies, provided that 
the displacement is now a measure of the degree of excitation of the appropriate vibrational 
mode of the sample. Acoustic resonances are then obtained when the elastic bodies 
vibrate in such a way that standing waves are set up (Some interesting papers dealing 
with diverse kinds of resonances in metals can be consulted in [42]). Since the profile of 
atomic phenomena involving a high number of quanta of excitation can be analyzed in 
the context of classical mechanics, cyclotron and electron spin resonances can be studied, 
in a first approach, in terms of the above model (see, e.g. the paper by A.S. Nowick 
in [42], pp 1-44, and references quoted therein). The quantum approach to the problem of 
spinning charged particles showing magnetic resonance is discussed in conventional books 
on quantum mechanics like the one of Cohen- Tannoudji et. al. [43]. 

We have also shown that the resonance phenomenon occurs in quantum decaying sys- 
tems. According with the Fock's approach, the corresponding law of decay depends only 
on the energy distribution function u(E) which is meromorphic and acquires the form of 
a FBW, bell-shaped curve. The introduction of e = E — iT/2, a complex eigenvalue of the 
energy, is then required in analyzing the resonances which, in turn, are identified with the 
decaying states of the system. The inverse of the lifetime is then in correspondence with 
T/2. In a simple model, the related exponential decreasing probability can be justified 
by introducing a complex potential U = Ur + Uj, the imaginary part of which is the con- 
stant — r/2. Thus, Uj works like a sink of probability waves. The situation resembles the 
absorption of electromagnetic waves by a medium with complex refractive index so that 
U = Ur + Uj is called an optical potential in nuclear physics. The treatment of complex 
energies in quantum mechanics includes non-square integrable Gamow-Siegert functions 
which are outside of the Hilbert spaces. In this sense, the complex-scaling method is use- 
ful to 'regularize' the problem by complex-rotating positions x and wavenumbers k. As a 
consequence, bounded and scattering states maintain without changes after the rotation 
while the Gamow-Siegert functions become square-integrable. Another important aspect 
of the method is that the positive real axis of the complex energy plane is clockwise rotated 
by an angle 29, so that resonant energies are exposed by the cuts of the corresponding 
Riemann surface. The method, however, also produces complex potentials. That is, the 
Gamow-Siegert functions are square-integrable solutions of a non-Hermitian Hamiltonian 
belonging to complex eigenvalues. Similar results are obtained by deforming the initial 
potential in terms of a Darboux transformation defined by a Gamow-Siegert function. In 
this sense, both of the above approaches could be applied in the studying of quantum 
resonances. A detailed analysis of the connection between the complex-scaling method 
and the Darboux-Gamow transformation is in progress. 



25 



Acknowledgments 



0R0 would like to thank to the organizers for the kind invitation to such interesting 
Summer School. Special thanks to Luis Manuel Montaho and Gabino Torres. SCyC 
thanks the members of Physics Department, Cinvestav, for kind hospitality The authors 
wish to thank Mauricio Carbajal for his interest in reading the manuscript. The support 
of CONACyT project 24233-50766-F and IPN grants CO FA A and SIP (IPN, Mexico) is 
acknowledged. 

References 

[1] Gamow G, Z Phys 51 (1928) 204-212 

[2] Breit G and Wigner EP, Phys Rev 49 (1936) 519-531 

[3] Bohm A, Gadella M and Mainland GB, Am J Phys 57 (1989) 1103-1108; 

Bohm A and Gadella M, Lecture Notes in Physics Vol 348 (New York NY: Springer, 
1981) 

[4] de la Madrid R and Gadella M, Am J Phys 70 (2002) 6262-638 

[5] de la Madrid R, AIP Conf Proc 885 (2007) 3-25 

[6] Civitarese O and Gadella M, Phys Rep 396 (2004) 41-113 

[7] Aguilar J and Combes JM, Commun Math Phys 22 (1971) 269-279; 
Balslev E and Combes JM, Commun Math Phys 22 (1971) 280-294 

[8] Simon B, Commun Math Phys 27 (1972) 1-9 

[9] Fock VA, Fundamentals of Quantum Mechanics, translated from the Russian by E 
Yankovsky (Moscow: URSS Publishers, 1976) 

[10] Hecht E and Zajac A, Optics (Massachusetts: Addison- Weslye 1974) 

[11] Jackson JD, Classical Electrodynamics, 3rd edition (USA: John Wiley and Sons, 
1999) 

[12] Landau LD and Lifshitz EM, The Classical Theory of Fields, 4th revised English 
edition, translated from the Russian by M. Hamermesh (Oxford: Pergamon Press, 
1987) 

[13] Brown JW and Churchill RV, Complex Variables and Applications, 7th edition (New 
York NY: McGraww Hill, 2003) 

[14] Fowles G.R., Introduction to Modern Optics, (New York: Dover Pub, 1975) 

[15] Siegert AJF, Phys. Rev. 56 (1939) 750-752 



26 



[16] Heitler W and Hu N, Nature 159 (1947) 776-777 

[17] de la Madrid R, Quantum Mechanics in Rigged Hilbert Space Language (Valladolid: 
PhD dissertation, Theoretical Physics Department, Universidad de Valladolid, 2003) 

[18] Dirac PAM, The principles of quantum mechanics (London: Oxford University 
Press, 1958) 

[19] Maurin K, Generalized eigenfunction expansions and unitary representations of 
topological groups (Warsav: Pan Stuvowe Wydawn, 1968) 

[20] Gelfand IM and Vilenkin NY Generalized functions Vol 4 (New York NY: Academic 
Press, 1968) 

[21] Fernandez- Garcia N and Rosas-Ortiz O, Ann. Phys 323 (2008) 1397-1414 
(larXiv:0810.5597i [math-ph]) 

[22] Bohm AR, Scurek R and Wikramasekara S, Rev Mex Fis 45 S2 (1999) 16-20 

[23] Fernandez- Garcia N, Estudio de Resonancias y Trans formaciones de Darboux- 
Gamow en Mecd-nica Cudntica (Mexico: PhD dissertation, Departamento de Ffsica, 
Cinvestav, 2008) 

[24] Nussenzveig HM, Causality and Dispersion Relations (New York: Academic Press, 
1972) 

[25] Fernandez-Garcia N and Rosas-Ortiz O, J. Phys. Conf. Series 128 (2008) 012044 
flarXiv: 0902 .40521 [quant-ph]) 

[26] Feshbach H, Porter CE and Weisskopf VF, Phys Rev 96 (1954) 448-464 

[27] Giraud BG and Kato K, Annals of Phys 308 (2003) 115-142 

[28] Sudarshan ECG, Chiu CB and Gorini V, Phys Rev D 18 (1978) 2914-2929 

[29] Mielnik B and Plebahski J, Ann. Inst. Henry Poincare XII (1970) 215 

[30] Berggren T, Phys Lett B 44 (1973) 23-25 

[31] Moiseyev N, Phys Rep 302 (1988) 211-293 

[32] Cannata F, Junker G and Trost J, Phys Lett A 246 (1998) 219-226; 

Andrianov AA, Ioffe MV, Cannata F and Dedonder JP, Int J Mod Phys A 14 (1999) 
2675-2688; 

Bagchi B, Mallik S and Quesne C, Int J Mod Phys A 16 (2001) 2859-2872 
[33] Fernandez DJ, Munoz R and Ramos A, Phys Lett A 308 (2003) 11-16 



(arXiv:quant-ph/02 12026) 



[34] Rosas-Ortiz O and Mu noz R, J Phys A: Math Gen 36 (2003) 8497-8506 
darXiv:quant-ph/0302190[ ) 



27 



[35] Rosas-Ortiz O, Rev Mex Fis 53 S2 (2007) 103-109 f!arXiv:0810. 2283 [quant-ph]) 

[36] Fernandez- Garcia, Rev Mex Fis 53 S4 (2007) 42-45 

[37] Munoz R, Phys Lett A, 345 (2005) 287-292; 

Samsonov BF and Pupasov AM, Phys Lett A 356 (2006) 210-214; 
Samsonov BF, Phys Lett A 358 (2006) 105-114 

[38] Mielnik B and Rosas-Ortiz O, J Phys A: Math Gen 37 (2004) 10007-10036 

[39] Fernandez DJ and R osu H, Rev Mex Fis 46 S2 (2000) 153-156 
( |arXiv:quant-ph/010209l"| ); 

Fernandez DJ and Rosu H, Phys Scr 64 (2001) 177-183 flarXiv:quant-ph/9910125) 

[40] Ramirez A and Miel nik B, Rev Mex Fis 49 S2 (2003) 130-133 
( |arXiv:quant-ph/0211048D 

[41] Sokolov AV, Andrianov AA and Cannata F, J Phys A: Math Gen 9 (2006) 10207- 
10227 

[42] Vogel FL (Editor), Resonance and Relaxation in Metals, Second (Revised) Edition, 
(New York: Plenum Press, 1964) 

[43] Cohen- Tannoudji C, Diu B and Laloe B, Quantum Mechanics, Vols 1 and 2, (New 
York: John Wiley and Sons, 1977) 



28 



